- (Brief) introduction to metabolomics
- Preprocessing of LC-MS data
- Normalization
- Annotation/identification
CSAMA 2019
Sample is dissolved in a fluid (mobile phase).
Mobile phase carries analytes through column (stationary phase).
Sample is dissolved in a fluid (mobile phase).
Mobile phase carries analytes through column (stationary phase).
Separation based on affinity for the column’s stationary phase.
Sample is dissolved in a fluid (mobile phase).
Mobile phase carries analytes through column (stationary phase).
Separation based on affinity for the column’s stationary phase.
HILIC (hyrophilic liquid interaction chromatography):
Hydrophilic, polar stationary phase.
Analytes solved in mobile phase.
Analytes separated by polarity: compounds with low polarity elute first, with high polarity later.
LC-MS data has 3 dimensions:
Analyze data along rt.
MSnbase: data import with readMSData.xcms: peak detection with findChromPeaks and algorithm-specific parameter object.cwp <- CentWaveParam(peakwidth = c(2, 10), snthresh = 5) data <- findChromPeaks(data, param = cwp) head(chromPeaks(data), n = 3)
## mz mzmin mzmax rt rtmin rtmax into intb ## CP001 114.0907 114.0899 114.0929 1.954 0.280 3.907 1559.829 1555.923 ## CP002 114.0913 114.0884 114.0929 5.860 4.465 8.650 1890.221 1885.757 ## CP003 114.0914 114.0899 114.0929 10.882 8.650 13.114 1950.953 1946.210 ## maxo sn sample ## CP001 584.9510 584 1 ## CP002 601.8881 601 1 ## CP003 691.9580 691 1
xcms: adjustRtime function with PeakDensityParam [Smith et al. Anal. chem. 2006] or ObiwarpParam [Prince et al. Anal. chem. 2006].xcms: groupChromPeaks with NearestPeaksParam [Katajamaa et al. Bioinformatics 2006] and PeakDensityParam [Smith et al. Anal. chem. 2006].xcms: groupChromPeaks with NearestPeaksParam [Katajamaa et al. Bioinformatics 2006] and PeakDensityParam [Smith et al. Anal. chem. 2006].
peak density approach (for a given m/z slice):
## DataFrame with 4 rows and 4 columns ## mzmed rtmed POOL_1 POOL_2 ## <numeric> <numeric> <numeric> <numeric> ## FT001 105.041814839707 167.96077773118 229.490739260736 3093.75184315684 ## FT002 105.041653033614 157.082078655729 4762.39872227773 6601.45091358641 ## FT003 105.069636149683 31.8109676438667 699.723986763237 1033.23232267732 ## FT004 105.11027064078 63.7514649484463 20211.2633706294 15839.5504368189
Annotate features based on m/z and:
MSnbase).xcms.